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1. Introduction 


At HERA, reactions of electrons or positrons and protons, ep —> eX, are probed at centre-of- 
mass energies of 320GeV. Two kinematic regimes are distinguished, depending on the negative 
momentum transfer squared Q 2 from the ingoing to the outgoing electron 1 . At high Q 2 > 4 GeV 2 , 
the process is referred to as deep-inelastic scattering (DIS), whereas at low Q 2 < 2 GeV 2 the pro¬ 
cess is called photoproduction (yp). In the analysis presented here [[I]], dijet production in both 
diffractive yp and diffractive DIS is studied. The reaction may be written as ep —> eXp, where 
X is a hadronic system which contains two jets. The outgoing proton is detected in the HI Very 
Forward Proton Spectrometer (VFPS) |^]|. 

QCD calculations for this process are based on 
diffractive parton densities (DPDFs). A corresponding 
factorisation theorem is known for the case of diffrac¬ 
tive DIS 0] but not for diffractive yp. The DPDFs de¬ 
scribe the probability to find a parton with longitudinal 
momentum fraction zip in the proton, given that there 
is a diffractive signature, characterised by a longitudi¬ 
nally momentum fraction 1 — xjp, momentum transfer 
t of the outgoing proton and a hard scale jl. They are 
folded with hard matrix elements to describe jet pro¬ 
duction at next-to-leading order in the strong coupling. 

A hard scale jd 2 = (Ej jet ) 2 + Q 2 is provided by the jet 
transverse momentum ( Ej jet ) and in the case of DIS 


by the momentum transfer Q 2 . 

Experimentally, the DPDFs are determined from 
inclusive diffractive DIS cross section measurements, 
where no requirements on the hadronic final state X are 
made. As suggested by the scheme shown in figure [I| 
the DPDFs are determined with the ad-hoc assumption 
that they factorise into a probability to find a colourless 

object IP in the proton and parton density functions, ascribed to the structure of IP. The probability 
or “flux factor” is taken to depend on xjp and t only, whereas the parton density functions of IP 
only depend on the variable zip and the hard scale /i. 

For the present analysis, the HI 2006 DPDF fit B [|j] is used to predict cross sections. The 
present measurements focuses on testing the factorisation assumptions both in the DIS and in the 
yp kinematic domain. Factorisation in diffractive DIS has been validated by experiment [|[ 
01 I®- In hadron-hadron collisions clear evidence for a suppression of the cross section as 
compared to the expectation by about one order of magnitude is observed [|lT[ |~i~2| |. For diffractive 
yp, dijet production has been measured at HERA [13, [M|, [l5[] . A suppression of about a factor 



Figure 1: Diffractive dijet production with 
a leading proton at HERA. 


of two is observed in the two independent analyses by the HI experiment, whereas in the ZEUS 
analysis no suppression is seen » 

1 Throughout this paper, the term electron or the variable e is used to denote both electrons and positrons, unless 

otherwise stated. 
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Figure 2: Control distributions of selected reconstructed quantities. 


2. Data Analysis 

The analysis [jl]] is based on events 
recorded with the HI detector 0] in the 
years 2006 and 2007, corresponding to an 
integrated luminosity of 30 pb 1 for yp and 
50pb -1 for DIS. DIS events are identified 
by the presence of an electron in the HI 
rear calorimeter (SpaCal) [[F7|], whereas yp 
events are identified by the absence of an 
electron in the SpaCal or the liquid-argon 
calorimeter. These requirements limit the 
accessible kinematic range to 4 < Q 2 < 

80GeV 2 for DIS and to Q 2 < 2GeV 2 for yp. 

The inelasticity y is also restricted, in order 
to ensure that both the scattered electron and the hadronic final state are well measured. Diffractive 
events with a leading proton are selected by requiring a tagged proton in the HI VFPS detector [^]. 
The VFPS has a high geometrical acceptance (full coverage in azimuth), but restricts the accessi¬ 
ble range in the forward proton longitudinal momentum xjp and in the momentum transfer at the 
proton vertex |f|. All selection criteria are summarised in table p]. The hadronic final state X is re¬ 
constructed from the detected tracks and calorimeter clusters, using an energy-flow algorithm. The 
four-vectors of the hadronic objects are boosted to the y*p rest-frame and jets are reconstructed 
using the inclusive kj jet algorithm [|T|] with Pj recombination scheme and distance parameter 
R = 1. Dijet events are accepted if there are at least two jets with transverse momenta of the 
leading (sub-leading) jet fulfilling the conditions Ep etl > 5.5GeV (Ej ,el2 > 4GeV). The jet pseu¬ 
dorapidity is restricted in the laboratory frame to — 1 < r] yf;n ’ 2 < 2.5, in order to ensure that the jets 



Photoproduction DIS 

Event 

Q 2 < 2GeV 2 4 < Q 2 < lOOGeV 2 

kinematics 

0.2 <y < 0.7 

Leading 

0.01 < xip < 0.024 

proton 

|f| <0.6GeV 2 

zip <0.8 

Dijets 

E* jen > 5.5GeV 

Ep etZ > 4GeV 
-1 < pj et1 ' 2 < 2.5 


Table 1: Analysis phase space. 
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are well contained 2 . The events are selected in an extended analysis phase space, in order to ensure 
that migrations near the phase space boundaries are be well controlled. Only after correcting for 
detector effects, the phase space is restricted to the boundaries described in the text. For unfolding 
from detector objects to the particle level, a regularised matrix unfolding technique is applied. The 
response matrix is constructed using the RAPGAP event generator interfaced to the GEANT-based 
simulation of the HI detector. Control distributions of reconstructed variables at detector level are 
shown in figure |j. The variables shown are: the local VFPS coordinates Xy PPS and Yy PPS , the 
proton fractional momentum loss xjp measured with the VFPS, the momentum transfer Q 2 (DIS 
events only), the leading jet transverse momentum E P et 1 and the average jet pseudorapidity (ij ' et ). 
All reconstructed variables are well described by the RAPGAP simulation. 

3. Cross section measurements 

Single-differential cross sections are measured as a func¬ 
tion of various variables. For yp it is particularly interesting 
to look at the variable x y , which is sensitive to the longitudi¬ 
nal momentum fraction of the photon entering the hard pro¬ 
cess. In a simplified view, the photon has a hadronic struc¬ 
ture at small x y but has no structure at x y = 1. As can be 
seen in figure [3, the shape in x y is well described by the 
NFO calculation, but the normalisation is not. This obser¬ 
vation is consistent with earlier measurements M 0- Hi- 
Integrated over the full phase-space, the cross section for di¬ 
jet production in diffractive DIS is measured to be 30.5 ± 

1.6(stat) ± 2.8(sysf)pb, in agreement with the NFO prediction 
of 28.3^g 1 4 4 (scale)+4Q(DPDF)±0.8(hadr). In yp, however, the 
cross section is measured to be 237 ± I4(stat) ±31 (syst) whereas 
the prediction is 430^gg 2 (scale)^gj(DPDF)± 13(hadr). The sys¬ 
tematic uncertainties, summarised in table stem in about equal 
part from the understanding of the VFPS detector, the hadronic 
energy scale, uncertainties of the RAPGAP simulation and the 
overall normalisation. 

Figure [| shows the cross sections in DIS and 
yp as a function of the variables zip, Ep etl and the 
difference in rapidity, |Atj|. In general, all single¬ 
differential cross sections are well described in 
DIS. In yp, the shapes are well described but the 
overall normalisation is off by approximately a 
factor of two. This effect also can be seen in fig¬ 
ure |5|, where the ratio of the measured to predicted 
cross section is depicted as a function of Q 2 . 

"The z axis is pointing along the proton flight direction. Polar angles 0 are measured with respect to the z axis. The 
pseudorapidity is defined as T] = — In tan (0/2). 


Systematic uncertainty 

DIS 

YP 

VFPS detector 

3.0% 

5.3% 

Hadronic energy scale 

4.4% 

7.2% 

RAPGAP model 

4.3% 

6.9% 

Normalisation 

6.0% 

6.0% 

Total 

9.1% 

12.8% 


Table 2: Summary of systematic uncertainties. 


^ HI VFPS data . AFG y-PDF 

” NLOH12006Fit-Bx0.83x(1+5 hadr ) 



Figure 3: Diffractive dijet pho¬ 
toproduction cross section as a 
function of the variable x y . 
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Figure 4: Diffractive dijet production with a leading proton; single-differential cross sections in DIS and in 
photoproduction. 


In the DIS regime, the ratio is compatible with one, whereas for photoproduction the ratio is close 
to one half. 


4. Cross section ratios 


In order to possibly reduce the impact of systematic uncer¬ 
tainties, ratios of cross sections in yp to DIS and double-ratios 
data to theory are also investigated. For the NLO prediction, 
the assumption is made that variations of the scale /r have to be 
applied consistently to the DIS and yp predictions. This leads 
to cancellations of scale uncertainties in the ratio. For the mea¬ 
surements, however, the systematic uncertainties do not reduce. 
The main reason are the model uncertainties, which are not cor¬ 
related between DIS and yp and thus are enlarged in the ratio. 
For the integrated cross section, the double-ratio of yp to DIS, 
data to theory is found to be 0.511 ±0.085(data) ' oo^ 2 (theory). 
This is consistent with an earlier HI analysis [ 13 ] of independent 
data. The present measurement, with a tagged forward proton, is 
free of proton dissociative contributions, as opposed to the earlier 
measurements. Hence one may conclude that proton dissociation 
certainly is not the dominant source of the suppression. 


^ HI VFPS data RAPGAP x 0.83 
S NLO H12006 Fit-B x 0.83 X (1+5 ) 



Q 2 [GeV 2 ] 


Figure 5: Ratio of measured to 
predicted diffractive dijet cross 
section as a function of Q 2 . 
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^ HI VFPS data RAPGAP 
NLO H12006 Fit-B x (1+S h ) 

3 H 2 =(E* ie, ) 2 +Q 2 — - Ji 2 =(E^ et, ) 2 +Q 2 /4 


^ HI VFPS data RAPGAP 
NLO HI 2006 Fit-B x (1 +5 hadr ) 

S GRVy-PDF -- AFGy-PDF 


^ HI VFPS data RAPGAP 

NLO HI2006 Fit-B x (1+5 hadr ) 

S GRVy-PDF . AFGy-PDF 


^ HI VFPS data RAPGAP 
NLO HI 2006 Fit-B x (1+5 had[ ) 

B GRVy-PDF -- AFGy-PDF 
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Figure 6: Ratio of yp to DIS diffractive dijet cross sections and double ratio of data to NLO prediction. 


Possible shape dependencies of the cross section ratios are also investigated. Figure |6| shows 
cross section ratios in the variables |Ar}|, y, zip and Ep et 1 . The double ratios are significantly 
different from unity, but no significant shape dependencies are observed. The largest deviation 
from a constant are observed in |A 771 , however the data are still compatible with a constant at a fit 
probability of 15%. 


5. Summary 


Diffractive dijet production with a leading proton is measured both in the regime of deep- 
inelastic scattering (DIS), Q 2 > 4GeV 2 , and in the photoproduction regime (yp), Q 2 < 2GeV 2 . 
The data exploit for the first time the HI Very Forward Proton Spectrometer, which has a large 
acceptance and thus small normalisation uncertainties. Single-differential cross sections, cross- 
section ratios and double ratios of yp to DIS are extracted. Next-to-leading order QCD (NLO) 
calculations are able to describe the measured diffractive DIS dijet cross sections, whereas they 
are off in normalisation in yp. The observed yp cross sections correspond to a suppression factor 
of 0.511 ± 0.085(data)+Q 0 0 -, 22 (theory), assuming that scale uncertainties of the NLO calculations 
largely cancel in the double ratio yp to DIS, data to NLO. This result is in agreement with earlier 
HI measurements [L3, @1- 
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